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Differential rotation on both components of the pre main-sequence 
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ABSTRACT 

We present the first measurements of surface differential rotation on a pre-main sequence 
binary system. Using intensity (Stokes I) and circularly polarised (Stokes V) timeseries spec- 
tra, taken over eleven nights at the Anglo- Australian Telescope (AAT), we incorporate a solar- 
like differential rotation law into the surface imaging process. We find that both components 
of the young, 18 Myr, HD 155555 (V824 Ara, G5IV + KOIV) binary system show significant 
differential rotation. The equator-pole laptimes as determined from the intensity spectra are 
80 days for the primary star and 163 days for the secondary. Similarly for the magnetic spec- 
tra we obtain equator-pole laptimes of 44 and 7 1 days respectively, showing that the shearing 
timescale of magnetic regions is approximately half that found for stellar spots. Both com- 
ponents are therefore found to have rates of differential rotation similar to those of the same 
spectral type main sequence single stars. The results for HD 155555 are therefore in contrast 
to those found in other, more evolved, binary systems where negligible or weak differential 
rotation has been discovered. We discuss two possible explanations for this; firstly that at the 
age of HD 155555 binary tidal forces have not yet had time to suppress differential rotation, 
secondly that the weak differential rotation previously observed on evolved binaries is a con- 
sequence of their large convection zone depths. We suggest that the latter is the more likely 
solution and show that both temperature and convection zone depth (from evolutionary mod- 
els) are good predictors of differential rotation strength. Finally, we also examine the possible 
consequences of the measured differential rotation on the interaction of binary star coronae. 
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1 INTRODUCTION 

Differential rotation is a key factor in the generation of stellar mag- 
netic fields. In the Sun it is differential rotation that transforms a 
large-scale polodial field into a stronger toroidal component. In the 
solar case this is thought to occur in the thin interface layer between 
the radiative core and the convective envelope called the tachocline. 
Recent Zeeman Doppler imaging (ZDI) studies of cool stars have 
found that this may not be the case for all stars. Regions of strong 
surface azimuthal magnetic field have been observed. The surface 
field of the Sun consists mainly of radial field and classical dy- 
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namo models suggest that any large-scale azimuthal field should be 
restricted to the tachocline. This has led to the suggestion that the 
dynamo in young rapidly rotating stars may not be restricted to the 
base of the convection zone a s in the Sun but be distributed through- 
out iDonati & Brownlll997r) . It is therefore important to measure 
stellar differential rotation to see if this provides any clues to ex- 
plain these observations. 

Techniques for measuring the strength of differential rotation 
by tracking the cool spots or magnetic regions have been developed 
from the results of Doppler imaging and Zeeman Doppler imaging 
(ZDI) studies. So far differential rotation has now been measured 
for more than 20 stars. These include stars of a ran ge of differen t 
surface temperatures and so convection zone depths. IB arnea J2005l) 



N.J. Dunstone et al. 



collates many of these results, finding that differential rotation rate 
is strongly related to surface temperature but only wealdy corre- 
lated with rotation rate. Near fully convective M-dwarfs (e.g. HK 
Aqr) exhibit very weak differential rotation, while early G-dwarfs 
have strong differential rotation. 

The stars considered in the lSamesI j2005l) study were all sin- 
gle young rapid rotators. Measurements of differential rotation have 
also been made for binary systems, again using Doppler imaging. 
The limited published sample currently contains the evol ved pri- 
mary staroftheRSCvnsy stems HR1099 dPetit et alj2004) and IM 
Peg JMarsden et al.l2007r) and the main-sequ ence secondary star o f 
the pre-cataclysmic variable star V471 Tau JHussain et al.ll2006h . 
All three stars are found to have weak differential rotation and in 
the case of V471 Tau it was consistent with solid body rotation. So 
the question remains whether it is the evolved nature of the stars 
or their membership of a binary system that is responsible for the 
observed decrease in differential rotation strength. In an attempt to 
disentangle these two effects we observed HD 155555 (V824 Ara), 
a pre-main sequence binary system consisting of two nearly equal 
mass components (G5IV -I- KOIV). 

In total eleven nights of observations were obtained of HD 
155555 (as detailed in j Q. A five night subset of the data was used 
in an earlier publication. iDunstone et al.l <2008r) (hereafter Paper I). 
This work focused on presenting the first magnetic maps of HD 
155555 and the longer term evolution of magnetic and brightness 
features. A binary Zeeman Doppler imaging co de was developed 
based upon the Doppler imaging code DoTS JCollier CameronI 
1 19970 and along the lines of previous binary Doppler ima ging 
codes (e.g. 



Vincent 



e lines or previous bmary Uoppler ima ging 
etalj 1 19931 , iHendrv & Mochnackil I2OO0I and 



IStrassmeier & RicduOOSl) . In Paper I we found that both compo- 



nents of HD 155555 have complex radial field magnetic topologies 
with mixed polarities at all latitudes. Also present were rings of az- 
imuthal field. HD 155555 therefore appears to share properties of 
both young single stars and the evolved primary stars of RS CVn 
binaries. 

In this paper we seek further clues as to the nature of dynamo 
processes occurring on this young binary system. We attempt to 
measure the differential rotation rates of both stars. In !j3]we use the 
subsequent five nights to produce an additional set of surface maps 
that are independent of those in Paper I. These two sets of indepen- 
dent maps are compared using the technique of cross-correlation 
in SJ4] We subsequently incorporate a latitudinal dependent shear 
into the binary imaging process in SJ5]and therefore use all avail- 
able data to measure the differential rotation strength. The results 
of our analysis are discussed in ^and we present our conclusions 
in 33 



2 OBSERVATIONS 

This paper is based upon spectropolarimetric observations made at 
the Anglo-Australian Telescope (AAT) using the University Col- 
lege London Echelle Spectrograph (UCLES) which was fibre fed 
by the SemelPol visiting polarimeter JSemel et al.ll 19930 mounted 
at the Cassegrain focus. Eleven uninterrupted nights of observa- 
tions were secured on our target star HD 155555 from 2007 March 
30 to April 09. Those observations obtained on the five nights 
March 31 - April 04 were analysed in Paper I. In this paper we 
combine all the observations to analyse the surface rotation prop- 
erties of HD 155555. The spectral extraction and reduction of the 
complete dataset was the same as that done for the subset of data 
used in Paper I and we refer the reader to this paper for a full ac- 



count of the exact process. In total we obtain a dataset of 106 Stokes 
I (intensity) and Stokes V (circularly polarised) spectra. 



3 COMPARING INDEPENDENT MAPS 

The 1.68 d orbital period of HD 155555 allows us to obtain full 
phase coverage in five nights. Given that we have a total timebase 
of eleven nights (March 30 - April 09), we can therefore create two 
entirely independent maps of both components of HD 155555. One 
set of brightness and magnetic maps has already been published in 
Paper I, corresponding to the 53 spectra taken between March 31 
and April 04 (hereafter referred to as dataset/epoch one). We there- 
fore use the subsequent five nights, April 05 to April 09 (hereafter 
referred to as dataset/epoch two), to create a second set of maps. 
Ten fewer spectra (43) are available for this second epoch which 
still provides excellent phase coverage. We note that the observa- 
tions taken on March 30 are not used in this part of our analysis. 
We prefer to have two equal five night datasets and therefore to ob- 
tain each complete map in as short a time as possible (to minimise 
the effect of differential rotation and evolution within each map). 
The effective time between the mid-points of the two epochs is five 
nights which is almost exactly three stellar rotations. 



3.1 Brightness maps 

As in Paper I, we u se the Doppler imaging code 'DoTS' 
JCollier Cameronll 19971) to map the surface brightness distribution 
of the surfaces of both stars using the Stokes I intensity spectra. The 
process is described in detail in Paper I and references therein. In 
Paper I standard tests were performed to explore the level of cross- 
talk in the maps of the two stellar components. We found this to be 
minimal due to the fact that for most of the orbit (approximately 
two thirds) the profiles of the two stars are cleanly separated in ve- 
locity space. The DoTS code has been extensively tested on binary 
systems with more entangled spectra than H D 155555, for exam - 
ple the contact system AE Phe was mapped bv lSarnes et al.l ( 120040 . 
The signal-to-noise (S/N) of this second five nights of spectra range 
between 95 - 210 in the peak order and are thus similar to that of 
the first. For brevity we do not show the obtained fits to the Stokes 
I data in this paper but instead refer the reader to Paper I where the 
fits to the first dataset can be examined. 

Observations are phased according to the ephemeris of 

IStrassmeier & Rice ( 2000) . see Table [T] Also listed in Table [T] are 
the system parameters for HD 155555 which were obtained in Pa - 
per I using the x^ minimisation technique of ISames et al.l ( 120000 . 
These orbital and physical parameters for the two stars are used in 
the Doppler imaging process. In Fig. [T] we present the brightness 
(spot) maps recovered from both epochs, for both stars. 

It is clear from Fig. [T] that both epochs produce very similar 
maps for the surface brightness distribution. The features described 
in Paper I are present also for the second dataset. On close examina- 
tion of the primary star map, a weak low latitude spot is recovered 
at phase <j) — 0.45 in the second dataset which has no counterpart 
in the first. The secondary star also has the major spot groups on 
the second epoch as were described in the first dataset. The two 
mid-latitude spots at phases — 0.7 — 0.8 are not as well defined 
on the second dataset as the first. In summary, the maps show little 
evidence of major evolution between the two datasets. The small 
differences that we do observe may be caused by one or more of 
the following factors: a) the presence of differential rotation; b) flux 
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Figure 1. Spots maps for both stars (primary - left, secondary - right) are plotted from the original Paper I epoch (top) and the new epoch (bottom). Vertical 
tick marks show the phases of observation. 



Table 1. Summary of the orbital and physical parameters for the HD 1 55555 
system taken from Paper I and references therein. 



Element 


Unit 


Value 


Orbital: 






P 


(days) 


1.6816463 


To 


(HJD) 


2446997.9102 


$0 




0.752474 


e 




0.0 


7 


(kms-i) 


3.72±0.02 


9 


m2/mi 


0.935±0.001 


Ki 


(kms-i) 


86.4±0.1 


K2 


(kms-i) 


94.7±0.2 


i 


(degs) 


50±5 (52 adopted) 


(i;sini)i 


(kms-i) 


34.9 


{v sin 1)2 


(kms-i) 


31.3 


Physical (assuming i = 52°): 




mi 


(Mq) 


1.054 


m2 


(Mq) 


0.986 


Rl 


(Rq) 


1.47 


R2 


(Rq) 


1.32 


TeB,l 


(K) 


5300±100 


TeS,2 


(K) 


505O±100 


log 51 




4.05±0.1 


log 52 




4.10±0.1 



emergence and diffusion; c) small changes caused by differences in 
the phase sampling of the two datasets. 



3.2 Magnetic maps 

We use our new ZDI code, 'ZDoTS', that was developed in Pa - 
per I (based upon the single-star ZDI code of lHussain et alfeOOOt) . 
to recover the magnetic maps of HD 155555 from the Stokes V 



spectra. ZDoTS is capable of modelling the contribution from each 
star to the combined Stokes V spectra and so makes use of the en- 
tire dat aset, including conjunction phases. As in the single star ZDI 
code of lHussain et al,l ( I200Q) . we model the local Stokes V assum- 
ing the intrinsic profile is Gaussian (which is fitted to the Stokes I 
profile) and constant over the stellar surface. We further assume that 
the weak field approximation holds. This is found to be the case for 
fields up to several kG from observati ons of very slowly rotating 
and w eakly magnetized Ap stars (see iDonati & Collier CameronI 
1 19971) , No assumption is made about the underlying photospheric 
temperature. The individual weights {w) assigned to each spectral 
line in the creation of the combined LSD profile is: w = gXd 
(where g is the magnetic Lande factor, A is the wavelength and d is 
the line depth). Inverse- variance weighting is then used to recover 
as optimally weighted and unbiased a representation o f the compos- 
ite Sto kes V profile as possible. W e refer the reader to lDonati et al,l 
1 1 19971) and lDonati & Brownl ( 11997.) for a full discussion of the im- 
plications of this approach. In addition to the tests of the ZDoTS 
code outlined in Paper I, further tests using synthetic data have been 
performed in iDunstona ( i2008i) . 

The magnetic polarity maps we recover yield only the mean 
magnetic flux density per surface resolution element, weighted by 
the mean photospheric brightness. Due to this fact the magnetic 
polarity signals from regions of low photospheric surface bright- 
ness (e,g, starspots) will be suppressed. The magnetic maps should 
therefore be considered as being representative mainly of the bright 
magnetic network. For the purposes of the present study this sup- 
pression has little influence on the differential rotation measure- 
ments because we are simply looking for identifiable small-scale 
magnetic structures that can act as tracers of the local differential 
rotation rate. In Fig, |2]we show the recovered radial field maps of 
HD 155555 for the two datasets. We do not show here the magnetic 
maps for the azimuthal field as the radial maps are sufficient to al- 
low the reader to appreciate the temporal evolution of the magnetic 
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Figure 2. Radial magnetic iield maps for both stars (primary - left, secondary - right) are plotted from the original Paper I epoch (top) and the new epoch 
(bottom). Vertical tick marks show the phases of observation. 



structures (the reader is again referred to Paper I if interested in 
examining the azimuthal field maps for the first epoch). 

As with the brightness maps, the radial magnetic maps of both 
components from the two epochs look very similar at first glance. 
On closer inspection, a number of features have changed their exact 
locations on the stellar surfaces, particularly on the primary star. 
As we shall determine in the following sections this is mainly due 
to the strong differential rotation for the magnetic features on the 
primary star. 



4 CROSS-CORRELATING MAPS 

The first technique that was employed to measure stellar differen- 
tial rotation using Doppler imaging involved cross-correlating im- 
age maps. This was done for the earl y differential rotation mea- 
surem ents on the KO dw arf AB Dor bv lDonati & Collier CameronI 
( 1 19970 and lDonati et alj al999) . The technique requires two inde- 
pendent maps that have been derived from observations spanning 
at least two rotations of the star. This is not difficult to achieve on 
the AB Dor as the 0.51 d rotational period allows almost complete 
phase coverage in a single night. Therefore even a two night observ- 
ing run could, in theory, be used for cross-correlation. However, 
in practise a couple of nights gap is normally left between epochs 
to allow enough time for measurable shear to occur. In compari- 
son, the 1.68 d rotational period of the components of HD 155555 
means that in 5 nights of observing we achieve complete phase cov- 
erage. From the last section we have two independent maps that are 

effectively separated by 5 nights (3 stella r ro tations). 

The cross-correlation JSimkinl 19741 and lTonrv & DavisI 19791) 
is performed on each latitude slice of the images. This is carried 
out in a way that accounts for the continuous nature of the stellar 
surface, i.e. that the surface maps in Figs. [T]&|2] wrap around. The 
resulting cross-correlation functions are shown in Fig. [3] and are 
examined for the characteristic signature of differential rotation. 



We find the peak of each latitude strip's cross-correlation function 
by fitting a parabola. This then allows us to obtain sub-pixel accu- 
racy in the resulting fits. Latitudes above 70° suffer from signifi- 
cant longitudinal smearing due to the poorer spatial resolution at 
high latitudes and so are not used. Similarly our maps contain no 
significant information below a latitude of -30° . The peaks of the 
cross-correlation functions are then themselves fitted with a solar- 
like sin"^ / differential rotation law (see below in ij5]( which is over- 
plotted on Fig. [3] 

The imaging code assumes that the rotational period of each 
star is equal to the binary orbital period. If this assumption of 
synchronous rotation is incorrect then the cross-correlation images 
should show an offset in phase. This would manifest itself as a ver- 
tical line offset from zero phase shift in Fig.[3l as an incorrect period 
would not be latitude dependant. From an examination of Fig. [3] it 
is obvious that all four maps show no evidence of uniform asyn- 
chronous rotation. The binary orbital period is the correct period 
at a particular stellar latitude. It is also clear that in all cases the 
equator is rotating faster than the poles. 

Out of the four cross-correlation images in Fig.[3]the only one 
that we are unable to fit with a differential rotation law was that 
derived from the primary star's brightness images. An examina- 
tion of the primary spot images in Fig. [T] shows why. Spots are 
present near the pole (above 60°) and at a low latitude band be- 
tween 5 and 35° but there are very few spot features at the in- 
termediate latitudes. The corresponding primary magnetic images 
(Fig. [2} however, have structure at all latitudes and so produce a far 
cleaner signal as shown in the bottom-left panel of Fig. [3] Here we 
find a differential rotation strength of ASl — 0.10 rad d~^, or an 
equator-pole lap time of 60 days. If we now compare this fit with 
the primary brightness map we can see that at the latitudes of 5- 
20°, within the low latitude band of spots, we find evidence of a 
similar phase shift. 

The cross-correlation of the secondary brightness map pro- 
duces a much better defined differential rotation signature than 
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Figure 3. The independent maps from the two epochs are cross-correlated to search for latitudinal dependant surface rotation. The two panels on the left are 
for the primary star, while the two on the right show the secondary star. The top plots come from cross-correlating the spot maps and the bottom plots from the 
geometric average of the cross-correlation images of the radial and azimuthal field maps. Plus symbols mark the maximum of the cross-correlation function 
for each latitude strip and the solid lines are the best fits of a simple solar-like differential rotation law for all but the primary spot cross-correlation image (see 
text). Note that the x-axis plotting phase has been reversed so as to present the reader with the more conventional diagram of differential rotation that would 
have been obtained if we had adopted stellar longitude. 



the primary star. Again this is because the secondary star has 
spots at a larger range of latitudes. It is well fitted by a sin'^ I 
law and results in a differential rotation measurement of Afl — 
0.04 radd^^, or a equator-pole lap time of 157 days. The cor- 
responding cross-correlation image using the secondary magnetic 
maps also produces a definite signature of differential rotation with 
Af7 — 0.06 rad d~^, or a near solar-like equator-pole lap time of 
105 days. 



5 IMAGE SHEAR 

Recent measurements of stellar differential rotation have been 
made using the sheared image technique, as developed by 
iDonati et alj feOOOh . Instead of cross-correlating two entirely inde- 
pendent maps, as in the last section, a latitudinal dependant shear 
is applied within the imaging code to model the effect of differen- 
tial rotation. As such each latitude strip has its own angular rotation 
rate given by: 



n{l) = fie 



Afisin I 



where (leg is the equatorial rotation rate and I is the stellar latitude . 
This was implemented into the DoTS code bv lSames et al.l J200lh 
to work for single stars. Here we implement this into the binary 



version of DoTS. This is achieved by first applying the latitudinal 
dependant shear to each vertex of our surface pixel grid for each 
star and then rotating the star in the co-ordinate frame of the bi- 
nary orbit. This change was also made to the our new ZDI code, 
ZDoTS, allowing us to obtain a separate measurement of the stellar 
differential rotation from the Stokes V spectra. 

The image shear method has many advantages over that of 
the original cross-correlation technique that we outlined in the last 
section. It works in data space directly and can therefore produce 
formal uncertainties on the differential rotation parameters. Many 
of the problems in fitting the differential rotation law to the peaks 
in the cross-correlation images also disappear. Furthermore, we do 
not require two complete and independent surface maps. Instead, 
all that is required is for at least one phase to be re-observed during 
a subsequent stellar rotation (preferably with both low and high lat- 
itude surface features). It is therefore possible to use quite incom- 
plete datasets (poor phase coverage) and will even work for stars 
with pathological periods near to a multiple of a day. Also the tech- 
nique is less sensitive to non-uniform latitudinal spot coverage (a 
problem we encountered with the brightness images of the primary 
star in f|4} as the shear is applied to the model rather than derived 
from the result. 

The values of ileq and Af7 are found by running a grid of 
models which sample combinations of the two variables. We solve 
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Figure 4. The results of the image shear technique for measuring differen- 
tial rotation using the spot (Stokes I) signatures. Differential rotation rate 
between the equator and pole is plotted as a function of absolute equato- 
rial rotation rate. The primary star is the upper plot and is found to have 
stronger differential rotation than the secondary star which is the lower plot. 
Contours showing the one parameter \-a, 2.6-cr (99%) and 4-cr (99.99%) 
confidence intervals are superimposed over the 2aS surface. 



for each set of rotation parameters for each of the two stars sepa- 
rately, as there is no reason to suspect that they should be signifi- 
cantly correlated. When obtaining the stellar parameters (displayed 
in Table [T] from Paper I) we set the code to aim for a very small 
value of x^ (in reality this is unobtainable) and then record the 
actual x^ achieved after a fixed number of maximum entropy itera- 
tions. The best value is then found by interpolating the two dimen- 
sional grid of parameters and searching for the global x^ minimum. 
We use a slightly different approach when obtaining the differential 
rotation parameters. 

The above technique does not take advantage of the im- 
age entropy information. The full relative posterior probability is 
exp(a6' — |x'^)' where S is the entropy of the final image and a 
is the Lagrange multiplier Therefore 2aS has the same units as x^ 
and when x^ is fixed can be used for purposes of determining error 
regions for the fitted parameters (Q,eq,AQ). So we choose an ob- 
tainable value of reduced x^ (e.g. x^=0-7 for the Stokes I data) and 
then record the value of 2aS once the code has converged. This 
is subtly different to the former technique and the probability con- 
tours are found to be smoother than that produced when using the 
minimum x^- This is probably due to the fact that the recovered im- 
age becomes dominated by spurious noise when pushing the image 
to small x^ and so is more susceptible to local minima. 

We apply the sheared image technique, as described above, to 
the entire 1 1 nights of observations (including the so far unused 
observations on March 30). In Fig. |4] we show the log posterior 
probability {2aS) surface defined by the sets of differential rotation 
parameters for each star from the intensity (Stokes I spectra) along 
with probability contours. As the two stars have different strengths 
of differential rotation it has been possible to display both stars on 
the same plot. Similarly in Fig. [5] we produce the same diagram 
but using the magnetic (Stokes V) spectra. The numerical results 
for the strength of the differential rotation (Af2) and the equatorial 
rotation rate (Qeq) are tabulated in Table |2] along with the results 
from fitting the cross-correlation functions in 21 The tilted ellipses 
of Figs.|4]c&|5]show that the parameters Af2 and Q^q are correlated. 
This is to be expected due to the fact that the shear rate is most 



Figure 5. The results of the image shear technique for measuring differen- 
tial rotation using the magnetic (Stokes V) signatures. Axis are as Fig. |4] 
Contours showing the one parameter l-tr, 2.6-(t (99%) confidence intervals 
are superimposed over the 2aS surface. The primary star is again found to 
have stronger differential rotation and is therefore the upper plot. 
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Figure 6. All four differential rotation measurements are plotted with l-cr 
uncertainties. Axis are as Figs. |4] and |5] The two squares are the measure- 
ments of differential rotation using spots and the two triangles are from 
magnetic features. Labelled lines join together measurements on each of 
the primary and secondary stars. The vertical dashed line shows the rota- 
tion rate corresponding to the binary orbital period. 



tightly constrained by mid-latitude features to which the Doppler 
imaging process is inherently most sensitive. 

The differential rotation result from the secondary star's in- 
tensity spectra of ASl = 0.0386 ± 0.0061 rad d^^ is very sim- 
ilar (consistent within the quoted uncertainty) to that produced 
from the cross-correlation technique of ASl — 0.04 rad d^^. This 
could have been predicted given that the sin^ law fit to the cross- 
correlation function of the secondary stars spot features, in the top- 
right panel of Fig. [3] was by far the most acceptable. The magnetic 
(Stokes V) spectra produce strengths of differential rotation con- 
sistently 30 % stronger than was found from the cross-correlation 
technique for both the primary and secondary components. 

We summarise the results of our analysis in Fig.[6]which com- 
bines all four of the differential rotation measurements from the 
sheared image technique. The primary star clearly has a stronger 
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Table 2. Surface rotation parameters for both components of HD 155555. Columns 2 and 7, labelled 'CCF', are the measurements from cross-correlating 
the independent maps produced from the two epochs in i|4] All other columns refer to measurements and associated uncertainties of the differential rotation 
strength (df2), the con'esponding equator-pole lap time ('Lap'), equatorial rotation rate (Hoq) and latitude of synchronous rotation (6^) produced from the 
sheared image technique in ^5\ 
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surface differential rotation rate from both the spot and magnetic 
features. Also the results obtained from the Stokes V spectra from 
both stars reveal that the differential rotation is approximately twice 
as strong from using the magnetic regions than the spots. 

Fig. [6] also shows that the equators of both stars rotate con- 
siderably faster than the binary orbital rate (shown as the vertical 
dashed line in Fig. [6j. This then raises the question of whether 
any stellar latitude is synchronously rotating with the orbit. In or- 
der to show this graphically we plot in Fig. [7] the differential rota- 
tion curves obtained for the two stars, for both the spots and mag- 
netic signatures. We clearly find that the poles are rotating sub- 
synchronously while the equators rotate super-synchronously in all 
cases. The exact synchronously rotating latitude for each differen- 
tial rotation measurement is given in Table[2] 

Curiously the primary and secondary stars share the same syn- 
chronously rotating latitude (47°) for the magnetic differential rota- 
tion. The resulting synchronously rotating latitude for the primary 
star's spot distribution is 39° and is thus the lowest of all four mea- 
surements. This may be interesting given that this latitude marks 
the upper boundary of the lower latitude band of spots seen on the 
primary star. 



6 DISCUSSION 

We have found that both components of HD 155555 have non- 
negligible surface differential rotation. This result was obtained 
from using the cross-correlation and the image shear techniques. 
There are two reasons why the cross-correlation technique (S|4j 
should probably be carried out whenever sufficient data of a star are 
available. Firstly, it provides a very literal interpretation of the data 
and so was useful in validating our implementation of the image 
shear technique to binary systems. Secondly, the cross-correlation 
technique does not assume a priori knowledge for the form of the 
differential rotation law which the image shear technique does. Be- 
yond this however, the image shear technique is superior and has 
many advantages which were outlined in ^and will be touched 
on again in i|6.4| We n o te tha t the direct spot tracking technique of 
ICollier Cameron et alj ( 120020 is another alternative which makes 
few assumptions but requires very good phase coverage and high 
signal-to-noise data. 



6.1 Internal velocity fields 

To summarise the results of ^ the primary star has stronger dif- 
ferential rotation than the secondary star and the magnetic regions 
produce stronge r signatures than do the spots. The latter trend was 
first revealed bv lPonati et Zl ( l2003h for AB Dor and HR 1099. The 



authors suggested that the spots and magnetic regions are anchored 
at different depths in the convection zones of these stars. The fact 
that the strength of the shear is so much greater in the magnetic re- 
gions suggests that the dynamos of these stars may not be active just 
at the interface layer between the radiative core and the convective 
envelope but instead be distribute d throughout t he con vection zone. 
Recent numerical simulations bv lBrown et al.l ( 120071) lend support 
to this idea by showing that global-scale toroidal and poloidal fields 
can be generated and maintained in the convection zones of rapidly 

rotat ing stars. 

iD onati et al.l ( 120030 attempted to quantify the nature of the in- 
ternal velocity fields in the rapidly rotating single stars AB Dor and 
LQ Hya by examining the temporal fluctuations in the differential 
rotation parameters (f2e<j,Af2). They found that: 

**eq ^ AZ\i2 ± ilsb 

where A depends on the assumed internal rotation model and on the 
internal stellar structure and fisi, is a constant, equal to the rotation 
rate the stellar convection zone would have if spinning as a solid 
body. 

For HD 155555 we only have a single epoch of differential 
rotation measurements however, given the binary nature of the sys- 
tem, we do have a very useful zero-point. We know that in the ab- 
sence of differential rotation, the rotation rate of each star would 
be equal to that of the orbital rotation rate. Therefore we can mea- 
sure A directly from Fig. |6] On this plot the gradient of the lines 
joining the zero point (where fl^q — Qsb and Afl — 0) with each 
of our differential rotation measurements will be equal to A~ ^. We 
obtain gradients in the range 1.6 to 2.5. lDonati et alj J2003b con- 
sidered two possible models for the angular velocity fields. Firstly, 
that angular rotation is constant with radius (like it is for the Sun) 
then our gradient (A^^) should be equal to five. Secondly that angu- 
lar velocity was constant along cylinders symmetric to the rotation 
axis, resulting in a gradient of approximately two (for stars such as 
AB Dor and the components of HD 155555). Therefore, like AB 
Dor, our results show that a model of internal rotation being con- 
stant along cylinders is most appropriate for the components of HD 
155555. Given their relatively rapid rotation this is what one would 
have expected. 

By ad opting this model of internal rotation, we can follow the 
example of iDonati et alj ( 120031) and calculate the range of possible 
latitudinal shears. At a given latitude the rotation rate will depend 
upon how deep the tracers of differential rotation are anchored in 
the stellar convection zone. The total depth of the convection zone 
of each star can be expressed as the fractional radius at the base 
of the convection zone, rc/R (where a star with rc/R=0 is fully 
convective and a star with rc/R=l has no convection zone). Then 
the maximum expected ratio between the shears we measure from 
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Figure 7. The differential rotation curves are plotted from the spots (left) and magnetic regions (right). The solid line is the primary star and the dashed line 
is the secondary star. The dot-dash vertical line corresponds to the orbital rotation period and thus shows that both stars are synchronously rotating with the 
binary orbit. 



the Stokes V and I data should be l/{rc/R)^. So if the observed 
value were to equal this maximum ratio then this would suggest 
that the magnetic regions are anchored near to the stellar surface 
while the spot features were anchored near the base of the convec- 
tion zone. The value of rc/R ca n be obtained using the evolution- 
ary models of ISiess et alj ( 1200(1) . The components of HD 155555 
have Vc/R = 0.68 for the primary and rc/R — 0.63 for the sec- 
ondary star. This corresponds to a maximum ratio of 2.16 for the 
primary star and 2.52 for the secondary star. Using our results for 
HD 155555, we find that the observed ratio between the shears from 
Stokes V and I is 1.8 for the primary star and 2.25 for the sec- 
ondary star We therefore find that HD 155555 is consistent with 
the adopted model of internal rotation. 



6.2 Comparison with other systems 

It is now possible to compare our results with those of both sin- 
gle rapidly rotating stars and binary systems. In the left panel of 
Fig. [8] we plot differential rotation rate versus effective tempera- 
tur e. A similar pl ot was first produced for single G, K and M dwarfs 
bv lBarnesI (|2005J). The authors discovered that for single stars the 

strength of differential rotation increases with increasing stellar sur- 

I ■ I ■ 

face temperature. More recently Reiners (2006) combined these re- 
sults with differential rotation measurements on F dwarfs derived 
from a Fourier an alysis of the ro tational broadening profile. An ex- 
trapolation of the lSarnesI ( 120050 power-law fit showed that the two 
techniques were in good agreement. The latest version of this dia- 
gram was made bv CameronI l|2008) who performed a re-analysis 
of the power law fit using all latest available differential rotation 
measurements. This is plotted as the solid line on Fig. [8] When 
we plot the differential rotation parameters for both components of 
HD 155555 we find that they effectively straddle this line. It there- 
fore appears that the strength of differential rotation on both stars 

is consistent with single stars of the same temperat ure. 

In additi o n to the single stars, collated in iBamej ( 120050 
and ICameronI (2008), we also consider the results from previ- 
ous Doppler imaging studies of binary star systems. Fig. [8] plots 
(as solid triangles) the differential rotation measurements of the 
evolved p rimary stars of th e R S CVn binaries HR 1099 and IM 
Peg (from lPetit et al.ll2004l and iMarsden et al.ll2007l respectively). 
The main-sequence component of the pre-cataclysmic variable star 



V471 Tau dHussain et al.l 120060 is also plotted as a filled cir- 
cle on Fig. [8] All three systems show weak differential rotation 
(An < 0.02 rad d~^). We note however that there is only a sin- 
gle epoch measurement for V471 Tau and this is consistent with 
zero {AQ. = 0.0016 ± 0.006 rad d"^) differential rotation. The 
recovered Doppler images suffer from a lack of spot coverage at all 
latitudes which may make them unreliable. We therefore encour- 
age future observations to confirm this measurement, but we do not 
consider V471 Tau further here. 



6.3 Tidal forces or internal structure? 

The RS CVn systems HR1099 and IM Peg have been repeatedly 
observed, so many measurements of differential rotation are avail- 
able. When compared to single main-sequence stars of similar mass 
HR 1099 and IM Peg have greatly reduced differential rotation. 
For example, HR 1099 has a mass of 1 Mq and yet exhibits a 
differential rotation rate approximately a third that of the Sun (or 
ap proximately an ei ghth that of the power-law fit shown on Fig. 
[8]l. IPetit et alj ( l2004h suggested that the reduced differential rota- 
tion was due to binary tidal forces. Over time these have acted to 
enforce synchronous rotation at all latitudes (i.e. to inhibit differen- 
tial rotation). The relatively strong differential rotation that we find 
on both components of the young HD 155555 system adds further 
complexity to the argument. We have shown that the mere fact that 
a star is a member of a close binary system does not necessarily re- 
sult in a reduced differential rotation rate. We are therefore left with 
two possibilities; that the timescale for tidal forces to suppress the 
differential rotation is long compared to the age of the HD 155555 
system, or that binarity is unimportant and evolutionary state is the 
dominant factor. We now examine the first of these possibilities. 

Our results for the young (18 Myr, Strassmeier & Rice 2000|) 
HD 155555 system show that tidal forces have indeed enforced 
synchronous rotation of both components, as shown in Fig.|7] How- 
ever, they have been unable to suppress the differential rotation, 
i.e. tides have not yet been able to synchronise rotation at all lati- 
tudes of the convective zone. Both components of HD 155555 are 
still contracting on to the main sequence. One may naively think 
that this would result in the stars rapidly spinning up and so might 
explain why the tidal forces are unable to effectively inhibit dif- 
ferential rotation. Given that HD 155555 is a close binary system 
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Figure 8. Left: Differential rotation is plotted as a function of stellar temperature. The solid line shows the power law fit of lCameroii (200s). Right: The 
X-axis now plots the fractional radius at the base of the convection zone (rc/R) as estimated from stellar evolutionary models (see text). Both plots include 
the same stars; filled symbols are members of binary systems, open symbols are single stars. The two triangles are the evolved RS CVn systems HR 1099 and 
IM Peg. The filled circle is the pre-cataclysmic system V471 Tau. Open squares are single G, K and M dwarfs, the 'star' symbols are F stars and the circled 
dot represents the Sun. The two large filled squares are the components of HD 155555. Error bars are plotted when the uncertainty significantly exceeds the 
symbol size. 



we should examine the efficiency of tidal forces. The synchronisa- 
tion timesca le {tsync) can be approximated as tsync—q^'^ {a/ RY' 
jZahnll 19771) . where q is the mass ratio. Using the parameters for 
HD 155555, we find a timescale of order tgync — 10^ yrs for both 
components. Given that this is very much shorter than evolution- 
ary timescales, it appears that tidal forces have ample opportunity 
to suppress differential rotation. Indeed, as the stars try to spin up, 
they will rapidly become resynchronised, with angular momentum 
being transferred from the stars' rotation to the binary orbit. The net 
result is that the stars will spin down rather than up. The fact that 
differential rotation is apparently uninhibited on the components 
of HD 155555 then places doubt on this mechanism also being re- 
sponsible for the low differential rotation observed on the evolved 
HR 1099 and IM Pe g systems. We futher note that theoretical work 
JScharlemannll 1982h predicts that tides should not inhibit differen- 
tial rotation on RS CVn systems. 

If binarity is unimportant then the strength of differential rota- 
tion is governed primarily by the stars evolutionary stage. HR 1099 
and IM Peg are both evolved giants and so have large radii (3.7 and 
13.3 Rq respectively) and correspondingly cool surface temper- 
atures. Therefore when we plot HR 1099 and IM Peg on the left 
panel of Fig. [8] (solid triangles) we find that they show strengths 
of differential rotation very close to (possibly slightly below but 
within the scatter) that of a main sequence single star with a sim- 
ilar surface temperature. This would seem to be a major clue as 
to the source of the low differential rotation. If tidal forces were 
solely responsible then there would be no reason to expect such an 
agreement with single stars. These giant stars have deep convection 
zones, like those of the cool K and M dwarfs, and therefore their 
internal structure may determine their differential rotation rates. 
In an attemp t to quantify this we use the evolutionary models of 
ISiess et alj ( 12000.) which quote the fractional radius at the base of 
the convection (rc/R). This constant was discussed earlier in %.U 

We use the evolutionary models to estimate the value of r^/R 
for all stars that were plotted on the temperature plot of Fig. [8] 
We then plot differential rotation as a fun ction of r^/R in th e right- 
hand panel of Fig.[8] Where the models of lSiess et al J ( 120001) did not 



follow the stellar evolution far enough the models of lClareg j2004) 
were used. The models are spaced in 0.1 Mq intervals and so we 
use the nearest one appropriate to each stars m ass which is obt ained 
from the literature. We only use F-stars from iReinersI ( 120061) that 
were defined as cluster members so that the cluster age could be 
used as a constraint to determine the stars position in the evolution- 
ary models. The single M dwarfs are particularly difficult to obtain 
values for Vc/R. This is because their evolutionary timescales are 
longer and so the depth of the convection zone continues to change 
for a longer time. Therefore we show approximate ranges for the 
value of rc/R where appropriate. 

In Fig. [8] we cover the entire range of convection zone depths. 
A strong trend of increasing differential rotation strength with de- 
creasing convection zone depth is found. This is very similar to 
the surface temperature plot which is not surprising as for main- 
sequence objects temperature essentially defines the depth of the 
convection zone. The position of the components of HD 155555 
therefore match well with the other single main- sequence stars. 
The RS CVn binaries HR 1099 and IM Peg move over to the left 
of the plot with values of rc/R = 0.16 and rc/R = 0.059 re- 
spectively. This illustrates their huge convection zone depths. At 
a specific temperature the convective depth is larger for an evolved 
star than for a main-sequence star. This would certainly seem to ex- 
plain the observed weak differential rotation. In fact HR 1099 and 
IM Peg have stronger differential rotation than would be expected 
for main sequence stars of similar fractional convection zone depth. 

We note that for consistency we have only considered the 
differential rotation results from stellar brightness maps (Doppler 
imaging). As found in !J5] the shear strength is often higher for 
magnetic regions than for star spots. Previous studies have shown 
that it is often easier to measure differential rotation using mag- 
netic maps because there is normally structure present at all stel- 
lar latitudes. It would therefore be interesting to extend the scope 
of this analysis an d consider results from ZDI also. For example, 
iMorin et al.l ( 12008!) found very weak differential rotation (Af2 — 
0.0063 ± 0.0004) on the fully conv ective M4 dwarf V374 Peg. 
In another recent studv lDonatil J2008l) found that the planetary host 
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star T Boo has very strong differential rotation {AQ, — 0.50±0. 12). 
This measurement is plotted on Fig.[8](solid diamond) and has a po- 
sition appropriate for its F7V spectral type. This is despite the fact 
that the stellar surface has apparently been f orced to rotate syn - 
chronously with the planetary orbital period l iCatala et al.l 120070 . 
So T Boo lends further support to the argument that tidal forces are 
unable to suppress differential rotation. 

6.4 Final maps 

The differential rotational parameters for both stars can now be in- 
corporated into the imaging process. We use the entire 11 night 
(106 spectra) dataset with the respective differential rotation pa- 
rameters for spots and magnetic features on each star to produce the 
brightness and radial magnetic maps shown in Figs. |9] and [TO] We 
note that the x^ achieved for these reconstructed images (x^=0.7 
for Stokes I and ;(^=1.0 for Stokes V) are essentially the same as the 
X^ for the individual images presented in ^ This would indicate 
that there has been little evolution in the spots or magnetic regions 
(other than the effect of differential rotation) over the eleven days. 
In order to show the effect of the differential rotation we also plot 
the primary radial field map obtained if we do not apply the differ- 
ential rotation law (top panel of Fig.fTOt. By including the effect of 
differential rotation the middle panel of Fig.llOlshows a higher level 
of detail. This is especially true at low latitudes where the effects 
of differential rotation are most noticeable. 

Looking back to Fig. [2] which shows magnetic maps produced 
from each five night subset of the data, we find that they contain 
more low latitude information than the top panel of Fig. [TO] (no 
differential rotation) but less than the middle panel Fig. 1101 (with 
differential rotation). This is what we would expect because the 
effect of differential rotation on only five nights is less than over 
eleven. Following the same argument we can now explain why the 
cross-correlation functions (shown in Fig.[3]( of the magnetic maps 
gave poor fits at low latitudes. Due to the relatively strong differ- 
ential rotation, even five nights is sufficient to smear the equator 
of the primary star by A</) = 0.11 and that of the secondary by 
Aff) — 0.07. Given that these are both larger than a whole reso- 
lution element (approximately A^ — 0.06) it is unsurprising that 
the cross-correlation function fails to follow the sin^ law at low lati- 
tudes. This once again illustrates the advantage of using the sheared 
image technique. 

Despite the improvement in the reconstruction of low latitude 
features, the maps presented in Fig. [lO] do not change the con- 
clusions of Paper I. It may be interesting to note however, that 
the small regions of flux, now well reconstructed between phases 
(j} = 0.3 — 0.5 and latitudes 20 and 30°, are in the same location 
as the gap in the ring of spots seen for primary star in Fig.|9] This 
apparent anti-correlation between the strengths of small magnetic 
regions and spot features may be an example of how the cooler spot 
temperatures suppress the contribution to the Stokes V spectra. 

6.5 Consequences for connecting field lines 

In Paper I we discussed the likelihood of field lines from one star 
connecting with the other. It was found from extrapolating the coro- 
nal field of each star separately (using the radial field maps) that the 
field complexity of the surface radial maps persists out to a signif- 
icant fraction of the binary separation (7.5 RqOT ~ 5.5 Rt, see 
Table [TJ and so the interaction between the two fields is likely to 
be complex also. Now that we know the strength of the differen- 
tial rotation on both stars, it is possible to speculate on the upper 





Figure 11. A cartoon sketch to show the effect of differential rotation on 
connecting field lines. The two components (primary - left, secondary - 
right) of HD 155555 and their separation are plotted to scale. A field line 
connects the equators of both stars. Initially at '1' the field footpoints are 
directly facing each other. At '2' we show the location of the footpoints 50 
days later when the combined angle swept out by the two footprints is 2it. 



limit of the lifetimes of such binary field lines and so the energy 
released from the rate of re-connection of the field. Let us con- 
sider a field line that is connected from a point on the equator of 
the facing hemisphere of one star to another on the equator of the 
facing hemisphere of the other star, a situation as sketched in Fig. 
nn Given that the two stars will be rotating in the same direction, 
but 180° out of phase, the two equatorial points will start to move 
away from each other due to differential rotation. The point on the 
primary star will sweep out a greater angle in the same time due to 
the stronger differential rotation (as shown in Fig.lUb. 

An interesting timescale often quoted for single stars is the 
equator-pole laptime. On this timescale a field line that connects 
the pole to the equator will have become wound around the star 
by a full circle (27r rads) and so stores significant magnetic energy. 
Similarly we can also ask; how long would it take the field points 
of our connecting field line to sweep out a combined angle of 27r 
rads? The answer is 50 days (illustrated in Fig.lllt which we note 
is comparable to the equator-pole lap-times of the primary star (44 
days) and the secondary star (71 days). Just like the equator-pole 
laptimes on a single star, the equator-equator field lines connecting 
the two stars gives us the maximum amount of stress on the field 
lines. We also note however that in the absence of any force con- 
fining such field lines to the equatorial plane the connecting loop 
will "balloon" up and over the stellar poles above a critical angle 
JLvnden-Bell & Boilvlll994h . This would further increase the rate 
of field reconnection events. 

The strong differential rotation provides extra stress on the 
field that would not be present on evolved binaries such as HR 1099 
and IM Peg with their weak differential rotation. We can there- 
fore speculate that the reconnection of these long binary field loops 
could significantly contribute to the X-ray luminosity of the HD 
155555 system and also to the frequency of large flares. The sce- 
nario outlined above caused by differential rotation is analogous to 
the mechanism thought to produce the so called 'superflares' on 
classical T Tauri stars. In this case field lines that connect the star 
with its surrounding disk become stret ched due to the fact tha t the 
disk is rotating slower than the star (e.g. lMontmerle et al.l200(]| . Fig. 
2). 



7 CONCLUSIONS 

We have measured the surface rotation properties of both stars in 
the HD 155555 system. A cross-correlation of the latitude strips of 
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Figure 9. Final brightness maps of the primary and secondary stars produced from all eleven nights and using the obtained values of differential rotation. 



both brightness and magnetic images, produced from two indepen- 
dent datasets, revealed the characteristic signature of differential 
rotation and confirmed that this can be modelled by a simple solar- 
like sin^ law. In order to better quantify the strengths of the differ- 
ential rotation we adapted the sheared image technique to binary 
systems. The rotational periods of both stars are found to be syn- 
chronously rotating with the binary orbital period. This confirms 
the fact that the system is tidally locked. 

Considering only the results from the brightness maps of the 
surface spot distributions, the primary has an equator-pole lap time 
of 80 d and the secondary star 162 d. These results show that 
the strength of differential rotation on both stars is not only non- 
negligible but are of similar magnitude to that found on single 
rapidly rotating stars of the same effective temperature. We used 
these findings to explore the relative effects of binarity and evo- 
lutionary stage on the strength of differential rotation. The large 
convection zone depths of the evolved giants in RS CVn systems 
is more likely to be the cause of their low differential rotation rates 
than the effects of tidal forces. 

Future measurements of differential rotation on the compo- 
nents HD 155555 are strongly encouraged. It would be very in- 
teresting to examine whether temporal changes in the differential 
rotation rate of the components of HD 155555 are correlated. This 
may indicate evidence of a joint activity cycle. 
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Figure 10. Final radial field magnedc maps. The top panel shows the primary radial field map obtained from all eleven nights of observations but without 
incorporating the measured differential rotation into the imaging process. This should then be compared with the middle panel where the differential rotation 
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rotation parameters. 
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